Extended Zee model for Neutrino Mass, Leptogenesis and Sterile Neutrino like Dark Matter 
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We propose an extension of the standard model with a L' ( 1 ) b -l global symmetry that accommodates radiative 
neutrino masses along with dark matter and leptogenesis. The observed matter antimatter asymmetry of the 
universe is generated through the leptogenesis route keeping the t/(l)B-L symmetry intact. The B — L global 
symmetry is then softly broken, providing the sub-eV neutrino masses. The model then incorporates a MeV 
scale sterile neutrino like dark matter. 
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I. INTRODUCTION 

In recent times, there have been growing interests on mod- 
els of neutrino mass, dark matter and baryon asymmetry yl 
|2[1. While the low energy neutrino oscillation data 1310,151] 
indicate that at least two of the physical left-handed (LH) neu- 
trinos have tiny masses and therefore they mix among them- 
selves, the galactic rotation curve |6], gravitational lensing [71] 
and large scale structure 18] strongly demand that there ex- 
ist non-baryonic dark matter in the present universe, which 
has only matter (visible plus dark) components. The antimat- 
ter components in the present universe are vanishingly small. 
This implies that the asymmetry between matter and antimat- 
ter components is maximal today. Currently this asymme- 
try has been precisely measured by the Wilkinson Microwave 
Anisotropy Probe (WMAP) ]|9J and is given by 
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where ub is the baryon density and riy is the primordial photon 
density. 

Despite the success of the standard model (SM), it can not 
explain any of the above phenomena: non-zero neutrino mass, 
existence of dark matter and matter-antimatter asymmetry of 
the universe. Explanation of all of these phenomena requires 
an extension of SM. The observed neutrino oscillation exper- 
iments deal only with the mass square differences of the neu- 
trinos, and hence, they leave an ambiguity on the nature of 
neutrinos to be either Dirac or Majorana. If the neutrinos 
are assumed to be Majorana particles, the sub-eV neutrino 
masses can naturally be generated through the celebrated see- 
saw mechanisms: type-1 1101 and type-II 1111 . In either case, 
the neutrino mass is suppressed by the scale of lepton (L) 
number violation. The smallness of the neutrino mass could 
also be explained naturally without invoking any large lepton 
number violating scale by the radiative mechanisms Ill2ll . The 
loop factor then introduces the suppression required to keep 
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the neutrino masses small. However, these models have the 
generic problem of explaining the baryon asymmetry of the 
universe, since the fast lepton number violation required to 
generate the neutrino masses, also erases any baryon asym- 
metry of the universe li 1 3il . In this note we aim to solve this 
problem and generate the neutrino masses and leptogenesis si- 
multaneously in a radiative model. The lepton asymmetry is 
generated without any B-L violation, similar to leptogenesis 
in models of Dirac neutrinos 11411 1511 . The present model also 
accommodates a sterile neutrino like dark matter candidate. In 
contrast to the generic radiative neutrino mass models, where 
the masses of the fields propagating in the loop are expected 
to be at the TeV scale, in the present model the masses of the 
fields propagating in the loop could be as high as lO'" GeV 
or so. As a result a leptogenesis could be possible from the 
out-of-equilibrium decay of these heavy charged scalars in the 
early universe. 

The paper is arranged as follows. In section II we propose 
a radiative neutrino mass model which has an exact B-L sym- 
metry to begin with. The B-L symmetry is then softly broken 
to give rise neutrino masses. Section III is devoted to estimate 
neutrino masses originating from the soft B-L violation. In 
section IV we calculate lepton asymmetry from a conserved 
B-L symmetry. In section V we estimate the lifetime of a de- 
caying sterile neutrino like dark matter. Finally section VI 
concludes. 



II. EXTENDED ZEE MODEL WITH CONSERVED B-L 
SYMMETRY 

Within the SM neutrinos are massless. This is because 
B — L is an exact symmetry of SM. So any radiative mech- 
anism within SM can not give rise to neutrino masses. One 
of the simplest extensions of SM is Zee model where SM is 
extended with a charged scalar, say r|^, and a second doublet 
Higgs, say (|)i, to generate neutrino masses through one loop 
radiative coiTection. The antisymmetric couplings of r|^ to 
Higgses i^^ah^'^a'^h) and leptons ifijV[^iiij) then together vi- 
olate lepton number by two units and neutrinos acquire masses 
through one loop radiative correction. Since the couplings 
of the Higgses and the charged leptons to the charged scalar 



is antisymmetric, the diagonal elements of the neutrino mass 
matrix vanishes identically in a flavor basis. The lepton num- 
ber violating couplings that appear in the neutrino masses, 
also give rise to lepton number violating interactions in the 
early universe. The observed neutrino masses would then im- 
ply that the lepton number violating interactions in the early 
universe is too fast, which will wash out fl3] any matter- 
antimatter asymmetry of the universe in the presence of the 
sphalerons before the electroweak phase transition. 

We propose an extension of the Zee model, which can 
simultaneously explain the observed neutrino masses, the 
baryon asymmetry of the universe and also provide a candi- 
date for dark matter of the universe. We extend the Zee model 
with an additional charged scalar %^ and three sterile neutrino 
like fields Nl, and impose a U{\)b-l global symmetry. The 
particle content of the model and their quantum numbers are 
presented in table |I1 The relevant terms in the Lagrangian are 



TABLE I: Particle content and their quantum numbers. 



Particle Content 


SU{3)cxSU{2)lxU{1)y 


U{1)b-l 


4 


(1,2,-1) 


-1 


^R 


(1,1,-2) 


-1 


(|>fl, h 


(1,2,1) 





X' 


(1,1,-2) 


-2 


na 


(1,1,-2) 





Nl 


(1,1,0) 


-1 



(^h can be given as 

((])„} = v„ and {(i?b) == v/, 
Then the two VEVs are related by 



(4) 



v2 + v^ = 174GeV, Va = vsm^ and Vi=vcosP, 



(5) 



where p = tan {I'a/vt,) 



III. SOFT B-L VIOLATION AND RADIATIVE NEUTRINO 
MASSES 

The [/(1)b-l symmetry is aUowed to be broken softly by 

Lsoft =m\'^x + h.c.. (6) 

As a result there is a mixing between r|^ and x^. In the flavor 
basis the mass matrix of r|^ and x is then given by 
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DiagonaUsing the above mass matrix we get the eigenvalues: 



^n'.x' 



{M\ +mI) ± J {Ml -M^f+4m^ 



(8) 



corresponding to the mass eigenstates r|' = cos 9r| + sin9x 
and x' = cos9x^ — sin9r|^, where 



then given as: 

^ 2 Mlx^x + Mlr\''^+fijx'iiLijL + {lJ)ab^Mb 

+hijn '^WiejR + Y^ti^dfaeiR + h.c. (2) 

where i,i — e,/j,z are family indices and a,b corresponds 
to two Higgs doublets. Since Nl is electrically neutral it 
can have Majorana masses Mat^. In Eq. ^ we have as- 
sumed that (^a couples only to leptons and (^b couples only to 
quarks (w\ apart from their self-interactions. However, this 
can be achieved by using an additional Z2 symmetry. We also 
assume that the couplings of leptons to (])„ is diagonal similar 
to ref. 1 16]. As a result there is no tree level flavor changing 
processes induced by Yij in the lepton sector. 

As the fields x and r|^ are charged they can not acquire 
any vacuum expectation value (VEV). The scalar potential in- 
volving the (|)„ and (^b then can be given as: 

+Xh\<^hf +Kb\<^a\^M^ +Kb\'^%f 

h.c] (3) 



-^[mbf- 



where M^jM^ > 0. The stability of the potential requires 
A-a, A-i > and Xab > —yKXb- The VEV of Higgses (|)„ and 
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Through the mixing between r|^ and x^ lepton number is 
violated by two units. As a result neutrino mass is generated 
through the one loop radiative correction diagram as shown in 
fig. [T] From fig. [T] the neutrino mass can be estimated to be 
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FIG. 1 : One loop radiative correction diagram for neutrino masses 
arises through the mixing between r| and %. 



{My),j = {Y[lMi - YfjMj) fij (m^iJabVb) ■ \I\ (10) 



where M, and Mj are the diagonal charged lepton mass matrix 
and the integral: 



IV. LEPTOGENESIS FROM A CONSERVED B-L 
SYMMETRY 



d q 



1 



1 



1 



1 



(271)4 ^2 - m2 ^2 _ m2 q2 _ j^l ^2 

'Ml\n{Ml/M^) +Ml\n{MyMl) +Ml\n{MllMl) 



16jt2 



(M2-M2)(M2-M2)(M2-M2) 



(1 



As the universe expands the temperature of the thermal 
bath falls. As a result r|^ will go out-of-thermal equilibrium 
below its mass scale. Note that r|^ has gauge interaction: 
r|^r|+ — > B^B^', where B^ is the U{1)y gauge field. How- 
ever, forMq ^ 10^" GeV the gauge interactions will remain in 
2 put-of- thermal equilibrium for a decoupled temperature, say. 



From Eq. ( fTOb it is clear that My is symmetric with respect 
to the family indices / and j and the diagonal elements of My 
are also zero. This is because of our assumption about the 
coupling of (|)a to the leptons. However, in general, this is not 
true. 

Ifwe assume that Mr|,M^ ^Ma,Mi,, then the neutrino mass 
is simply 
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and Fij = fij?>ij, where 



';' 



cotp 



(12) 



ln(M2/Mg) 

M2 



ln{Ml/Ml) 



Ml 



(13) 



To — Mf^/lO. The partial decay width of r|" 
given as 
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where the family index /, j = e,/J,T. At a cosmic temperature 
T ~ Mr|, if Fri fails to compete with the Hubble expansion 
parameter 
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where §* = 106 is the number of relativistic degrees of free- 
dom, then r| goes out-of-thermal equilibrium. From equations 
(fTTT i and ( fTST i then we find that the out-of-equilibrium decay 
of T]^ occurs for 



Mr\ > 2.77 X 10'°GeV (t^ 



(19) 



Note that decay of r] produce a pair of lepton and anti- 



Mj - M] 
^" " M} 



(14) 



We now estimate the magnitude of front factor Mq by using 
the following sample set: 

Mr, = 3 X lO'^GeV , M^ = 5 x 10"^GeV and M« = SOOGeV 

(15) 
We also choose the mas dimension coupling to be /j = 10^^ 
GeV and the soft B-L violation scale m = 10^ GeV. Then we 
get the symmetric neutrino mass matrix 



(Mv),-, =0.3 eV cotp /;-,5,- 



(16) 



The coupling constants Fij — fijdij can be appropriately cho- 
sen to explain the current neutrino oscillation data. Similar to 
the other neutrino mass models the coupling constants Fjj are 
constrained by the flavor changing processes like /j ^ e + y, 
T -^ fi + j, etc.... 

Note that in contrast to the original Zee model here the 
charged scalars flowing through the loop are super heavy. 
Therefore, they remain in out-of-thermal equilibrium for a 
while in the early universe, and hence, can generate a net lep- 
ton asymmetry consistently as discussed below. 





FIG. 2: Tree level and one loop self energy correction diagrams 
whose interference generates a net CP asymmetry. 



lepton. Therefore, the decay of r| does not produce any lepton 
asymmetry. However, if there are at least two r| ^ fields then 
there can be CP violation in the decay if r| fields. In the pres- 
ence of their interactions, the diagonal mass Mi in equation 
(|2]i is replaced by 



^U^l)cd^d + i^:■)H^-)cd^*d 
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(21) 



r;..^ii.' 



and gcc = Tt^ccMt], 



with Fr|(c = Fr|^ . Similarly, the interaction term /Jab'^'^a'^h in 



equation ^ should be replaced by {/Jab)c^c<^a'^b = l^c^c<^a'^b 
and {^^ab)d^d'^a'^b = /^d'^d'^a'^b^ where we have suppressed the 
symbols "ab" in /j. Now the absorptive part of one loop self 
energy diagram for rj^ -^ rj^/ can be given by 



Tn., A^% = ^ PcP*, + M^^M^, Y. hnjh*sj • (22) 



Diagonalising the mass matrix (|2T] ) one will get two mass 
eigenvalues corresponding to the two mass eigenstates \(/j and 
\(/2 . Note that the mass eigenstates \|/[^ and \(/[" (similarly xi/j" 
and "^2) ^re not CP conjugate of each other even though they 
are degenerate mass eigenstates, while r|j^ and r|j" (similarly 
r|+ and rjj^) are CP conjugates of each other Therefore, the 
decay of the lightest \(/^, say i|/]^ can generate a net CP asym- 
metry through the interference of tree level and self energy 
correction dia gram as shown in fig. (|2]). The CP asymmetry is 
then given by lUTT . 
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(23) 



Note that there is no radiative correction diagram. The decay 
of vi/]*^ does not violate L-number, since the decay of r|* does 
not violate lepton number. Therefore, the out-of-equilibrium- 
decay of \\if does not produce any L-asymmetry. However, 
the decay of \(/]^, below its mass scale, generates an equal B-L 
asymmetry between Nl and en due to the CP violation i^. 
The B-L asymmetry stored in e^ is then transferred to ei 
through the f-channel process e^e^ ^^ '^'^1 ^^ ^l^l ^^ shown 
in the figure ^. These interactions will remain in thermal 
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FIG. 3: The L-number conserving process which transfer the B-L 
asymmetry from right handed sector to the left-handed sector. 



equilibrium for all three generations of charged leptons below 
10^ GeV and hence there will be an equal amount of eg and 
ei asymmetry. The B-L asymmetry in e^ will be converted to 
the baryon asymmetry of the universe before the electroweak 
phase transition when the sphaleron processes are in thermal 
equilibrium, while an equal asymmetry will remain in Nl- 
The two asymmetries will equilibrate with each other after the 
electroweak phase transition when the sphaleron processes go 
out-of-thermal equilibrium. 

The final baryon asymmetry thus generated can be given as 
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where K = T\/H, Y\ being the decay rate of \|/i, measures 
the deviation from equilibrium. If K <^ \ then the baryon 
asymmetry is simply {ns/s) ^ £i/g*- On the other hand, if 
^ > 1, then the final baryon asymmetry suffers from a sup- 
pression of l/K{\nK)'^-^. If we assume that the CP violation 
is maximal, then by substituting M^J{M^^ ~^r\i) — '^i^) 
and /J2/a'i — 0(1) in equation (|23l ) we get a CP asymmetry 
El ^ 10^^ for huj ~ h2ij ^ 10^^. Smaller values of Ei can 
be conspired if we assume non-maximal CP violation. If we 
further assume that K —I. Ox 10^, then the suppression factor 
will be 3 X 10^^. As a result we will get a net baryon asym- 
metry hb/s ~ 10^'". The observed baryon asymmetry is then 
given by r| = ns/riy — l{nB/s). An exact value of the suppres- 
sion factor can be obtained by solving the required Boltzmann 
equation numerically which is beyond the scope of this paper 



A. Washout Constraints 

Below the mass scale of lightest r|^, the interaction 
h*:(f\^'ejRNiL is already gone out-of-thermal equilibrium. 
Therefore, there is no direct transfer of B-L asymmetry stored 
in Nl to e^. However, in a thermal bath the scattering: 
NiL^jR -^ Vjieji can occur through the mixing between r|^ 
and %^ as shown in fig. (|4|i. Note that this process violate 
lepton number by two units and therefore it is suppressed by 
the large mass scale of r|^ and x^. 
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FIG. 4: Scattering of Ni and e« through the mixing between r\ and 



To ensure that, we compute the scattering cross-section at a 
temperature below the mass scale of lightest r|^: 



(a|v|) = 
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where we have assumed that Ee — En = E. Thus from the 
above equation we see that below the mass scale of lightest r|^ 
and x^, the cross-section goes as: (a|v|) °^ ni^/M^, assuming 
that Mt| « M^ = M. The corresponding scattering rate for Nti, 
i ^ e,/j,z can be given as: 



^NL^n,^ {a\v\) 



(26) 



where n^ is the equilibrium number density of Nl and it is 
given by 
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The out-of-equilibrium of the scattering process: NilSjr 
ViL^jL then requires 



M{T:^M), 



(28) 



where H is the Hubble expansion parameter From equation 
I we get a constraint on the soft B-L violation scale to be 
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where we have used g^: — 106. Thus from the above equa- 
tion we get a constraint: m < 0(10'") GeV in order that the 
scattering will remain in out-of-thermal equilibrium above the 
electroweak phase transition. 



B. Low reheat temperature and Viability of leptogenesis 

It is believed that the universe has gone through a period 
of inflation and then reheated to a uniform temperature TJeh- 
If the corresponding theory of matter is supersymmetric then 
Tieh is highly constrained by the success of Big-Bang nucle- 
osynthesis, which could be spoiled by the overproduction of 
gravitino during the radiation dominated epoch |18]. For a 
gravitino mass of 0(100 GeV -1 TeV), a conservative up- 
per bound on Tj-^h reads 10^^^ GeV Ill9ll . If r^h is the maxi- 
mum temperature during reheating then it is difficult to create 
sufficiently high number densities of GUT gauge and Higgs 
bosons including the super heavy r|^ and x^. However, as 
discussed in refs. ll20l ulll . after the inflationary era the tem- 
perature does not rise instantaneously to Treh, but rises ini- 
tially to a maximum temperature Tmax and then falls to Tj-eh. 
It is argued that Tmax can be as high as lO^Tjeh 12 111 . As a re- 
sult the super-heavy charged particles r|^ and x* can be eas- 
ily produced through the gauge interactions: B^B^^ — *■ r|+r|^ 
and B^iB^^ -^ X^X ■ Even though T^nax is quite higher than 
Tieh, the gravitinos are mostly produced at T^^h ll22ll and there- 
fore, Tmax 3> Ti-eh is no more harmful for the success of Big- 
Bang nucleosynthesis. Subsequently the CP violating decay 
of these charged particles can produce lepton asymmetry con- 
sistently as discussed above. 



matter It will decay through the three body process: A^i -^ 
e^e+V, where the neutrino can be in any generation. Since 
this process violate B-L by two unit, it is naturally suppressed 
by the large mass scale of % and r|. The decay rate can be 
given as: 
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where hie is the coupling of A^i to e+ and r| , while fg-c is the 
coupling of x^ to e^ and Vx. Taking the B-L violating scale 



m = 10^ GeV, Mr, = 3 x 10'" GeV and My,^5x 10"' GeV, 
as taken previously, the life time of A^i is found to be 



Zn, = 0.88 X lO^^'Sec 



20c 



0(1) 
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Thus we see that Xni ^ To, where To = 0.18 x lO'^ Sec, is 
the age of the universe. Thus A^i can be a candidate of dark 
matter. 



VI. CONCLUSIONS 

We proposed an extension of the Zee model with a con- 
served B-L global symmetry. The B-L symmetry is softly 
broken to give rise to the neutrino masses through one loop ra- 
diative correction. In contrast to the Zee model, in the present 
case the charged scalars flowing in the neutrino loop are su- 
per heavy. Therefore, these charged scalars could depart from 
thermal equilibrium in the early universe. As a result the CP 
violating decay of the super heavy charged particles, namely 
r|-, could generate a net baryon asymmetry through the lepto- 
genesis route. Recall that the lepton asymmetry is generated 
without any B-L violation. This model then accommodate a 
sterile neutrino like dark matter A^i as its three body decay 
A^i -^ e^e+v is suppressed by the large mass scale of r| and %. 
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